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I. INTRODUCTION 

Weakly Interacting Massive Particles (WIMPs) are among the most popular candidates for 
the relic cold dark matter (DM). There is some revival onirterest_iii_tlie_WIMP-nucleus spin- 



dependent interaction 



[see e.g. 



Tom bot 



1 theoreticalfsee e 



0,0, ll,|l2, 13, 14, 15,0, 12, 18, Q 



) and experimental 



) points of view. There are some proposals 



aimed at direct DM detection with relatively low-mass isotope targets jS 
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as 



well as some att emp ts to design and construct a DM detector which is sensitive to the nuclear 



recoil direction 



; emp ts t 
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26|. Low-mass targets make preference for the low- 



mass WIMPs and are more sensitive to the spin- dependent WIMP-nucleus interaction as well 
Q, 3, y, 0,0, HQ- On the other hand, WIMPs with masses about 100 GeV/c^ follow from 
the results of the DArk MAtter (DAMA) experiment. This collaboration claimed observation of 
;he first evidence for the dark matter signal due to registration of the annual modulation effect 
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2^. 



Aimed for more than one decade at the DM particle direct detection, the DAMA 



experiment with 100-kg highly radio-pure Nal(Tl) scintillator detectors successfully operated 
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till July 2002 at the Gran Sasso National Laboratory of the I. N.F.N. On the basis of the results 
obtained for over 7 annual cycles (107731 kg-day total exposure) the effectiveness of the WIMP 
model-independent annual modulation signature was demonstrated and the WIMP presence in 
the galactic halo is strongly supported at 6.3 a C.L. |3l| . 

The goal of this review (being a continuation of our previous review paper Q) is to complete 
our nuclear physics consideration of the spin-dependent (or axial-vector) interaction of dark 
matter particles with nuclei. This type of interaction of the DM particles is important for the 
following reasons: (i) the spin-dependent interaction of the DM particles provides us with twice 
stronger constraints on the SUSY parameter space in comparison with the spin-independent 
interaction; (ii) in the case of spin- dependent interaction of heavy WIMPs with heavy target 
nuclei the so-called long g-tail behavior of the relevant form-factor allows detection of large 
nuclear recoil energy due to some nuclear structure effects; (iii) it is worthwhile to note that 
by relying only upon the scalar interaction of the DM particles, which seems to be strongly 
suppressed, one might miss a DM signal 0|. However, by a simultaneous study of both spin- 
dependent and spin-independent interactions of the DM particles with nuclei the chance for 
observing the DM signal is significantly increased. 

A dark matter event is elastic scattering of a relic WIMP (or neutralino) x with mass 
m-^ from a target nucleus A producing a nuclear recoil E-^ which can be detected by a suitable 
detector. The differential event rate in respect to the recoil energy is the subject of experimental 
measurements. The rate depends on the distribution of the relic WIMPs in the solar vicinity 
f{v) and the cross section of WIMP-nucleus elastic scattering 
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The differential event rate per unit mass of the target material has the form 



dR ^^'max da 



A 



I dvf{v)v-^{v,q 



2^ 



dE^ Jv^in dq 

We assume WIMPs (neutralinos) to be a dominant component of the DM halo of our galaxy with 
a density = 0.3 GeV/cm^ in the solar vicinity. The nuclear recoil energy En = q^/{2MA) 
is typically about 10~^m^ and Nt is the number density of a target nuclei with mass Ma- 
The WIMP-nucleus differential elastic scattering cross section for spin- non- zero (J 7^ Ojnuclei 



contains coherent (spin-independent, or SI) and axial (spin-dependent, or SD) terms ll. l39l l40l|: 
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Here ua = — — is the reduced WIMP-nucleus mass, o"4(0) = —lA^Cn] is the spin- 
rUy^ + MA ' siv y ^ L OJ 

independent WIMP-nucleus total elastic cross section at = and -Fgj(g^) is the normal- 
ized (Fgi(O) = 1) nozero- momentum-transfer spin-independent nuclear form-factor. The g^- 
dependence of the SD cross section is governed by the spin-dependent structure function S^Yiil) 
of Engel et al. i Q 

SiM=S''{q) = E i\{m\rf'{q)\\N)\' + \{N\\Cl{q)\\N)\'). (3) 

L odd 

The transverse electric T^^^{q) and longitudinal C^{q) multipole projections of the axial vector 
current operator are given in the Appendix. 

Relic WIMPs in the halo of our Galaxy have a mean velocity of (v) ~ 300 km/s = lO^^c. 
When R <^ l/^max (^^ Q'max-R ^ 1) whcrc R is the nuclear size and gmax = 2/iyif is the maximum 
of the momentum transfer in the ^-nucleus scattering, the spin structure function S^{q) reduces 
to (so-called zero transfer momentum limit) 

^^(0) = ^l(^llE^(«o + air>.||A)p = ^J(J + l)A^ with a = ^ + 

i 

Here (Sp,„) is the proton (or neutron) spin averaged over nucleus A, an and ap are the effective 
spin WIMP-neutron and WIMP-proton couplings that contain details of the supersymmetric 
model, as well as the quark spin content of the proton and neutron. 

As m-^ increases, ~ 1/g (the product qR starts to become non-negligible) and the finite 
momentum transfer limit must be considered for heavier nuclei. The formalism is a straight 
forward extension of that developed for the study of weak and electromagnetic semi-leptonic 
interactions in nuclei. Here we follow the definitions of 0,^^. With the isoscalar spin coupling 
constant = ctn + ap and the isovector spin coupling constant ai = ap — a„ one can split the 
nuclear structure function S'^^q) into a pure isoscalar term, S^Q^q), a pure isovector term, S^i{q), 
and an interference term, S^i^q), in the following way: 

S^'iq) = alS^,{q) + afS^M + aoa^M- (4) 

These three partial structure functions contain expectation values of operators of the form 
ji(gr)[Y^cr]'^^^, which depend on spatial coordinates and nucleon spins. The relations 



^ol,(o) = ci-msp) + (s„))^ stM = cimsp) - (s„))^ (s) 

S^,{0) = 2C{J){{Sl) - {SD) with C{J)- 



4n J 
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connect the nuclear spin structure function S^{q) at g = with proton (Sp) and neutron (S„) 
spin contributions averaged over the nucleus. In relations (0) the normalization coefficient 
C(J) > and therefore 5'oo(0) > and 5'^i(0) > 0. These three partial structure functions 
SfAq) allow calculation of spin-dependent cross sections for any heavy Majorana particle as 



42|. 



well as for the neutralino with arbitrary composition 

In this paper we extend our consideration of the modern nuclear spin structure calculations 
involved into the problem of the direct dark matter search. The calculations of the proton 
and neutron spins (Sp(„)) averaged over all nucleons in the nucleus A, which are relevant to the 
zero-momentum neutralino- nucleon spin cross sections, are considered in our previous review [3]. 
The cross sections at zero momentum transfer show strong dependence on the nuclear structure 
of the ground state j^]. Here we discuss the calculations of the spin structure functions in the 
finite momentum transfer approximation. We also touch upon the level of accuracy of these 
calculations. Finally, (in the last section) we briefly discuss a new approach to data analysis 
in the finite momentum transfer approximation directly in terms of the effective spin nucleon 
couplings ao,i together with the scalar WIMP-proton cross section cTgj(O) at g = 0. Contrary to 
other possibilities, (see for example, ^^i^ procedure is direct and relies as much as possible 
on the results of the most accurate nuclear spin structure calculations. 



II. NON-ZERO MOMENTUM TRANSFER LIMIT 

To the best of our knowledge, the finite, non-zero momentum transfer calculations of the 
spin nuclear structure functions S^{q) have been performed for the set of isotopes given in 
Table m The zero- momentum transfer limit (mostly quenching) is also investigated for Cd, Cs, 
47,Q,0, for hydrogen, ^H, H Q> helium, ^He, l^j, chlorine, ^SQ, Q and 



Ba and La 
calcium, "^^Ca, [3c 

General discussion of the nuclear model approaches to calculation of the spin characteristics 
. ike spins averaged over the nucleus proton (neutron) (Sp(„)), one can find in our previous paper 
~\. In this paper all available sets of the spin structure functions are given either in the form 
of explicit functions or as useful analytical parameterizations of the accurate numerical results, 
or only graphically (as pictures from original papers). 

As already noted in the introduction for quite heavy WIMPs and sufficiently heavy nuclei. 
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TABLE I: List of isotopes with available spin structure functions, S^{q), at g > 0. 



A 




/T-ULliUio ciliU. i t!lt!i tJllL-ct o 1 


19 


Fluorine, "'^^F 


Vergados et al. [26, 28, 44] 


zo 


oOQluni, rsu, 




97 
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9Q 


OlilLUli, Ol 


XVcoocli t;L di. 
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"Pirl" Q Gcm TTi 

± UldoolU-lli, Xv 


F,na-pl pf al [49] 
XL/li^cl clj di. 


73 


Germanium, ^^Ge 


Ressell et al. [40] 
Demitrov et al. f45l 


aO 


IN iUUlUlli, IN U. 


TTnrrpl pf al Mfi! 




' 1 Vul 1 n T*m Tn 123 ' | W 
X cll U.i 1 U.lil; -LC 
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125 


Tellurium, ^^^Te 


Ressell and Dean f4l1 


127 


Iodide, 


Ressell and Dean [41] 


129 


Xenon, ^^^Xe 


Ressell and Dean [41] 


131 


Xenon, ^^^Xe 


Engel [1] 

Ressell and Dean f4l1 

Nikolaev and Klapdor-Kleingrothaus [4^] 


207 


Lead, 207pb 


Vergados and Kosmos [44. 481 



the dependence of the nuclear matrix elements on the momentum transfer cannot be ignored 
even if the WIMP has energies as low as 100 keV. For example, if ^ m^, the reduced mass 

almost reaches (/x^ — '^a)- It is rather popular (and simplest) assumption that the 
WIMPs have a Maxwellian velocity distribution in the halo of our Galaxy. Some WIMPs will 
possess velocities significantly greater than (v) ~ 10~'^c. A maximum velocity of fmax — 700 
km/s (slightly greater than the Galactic escape velocity and more than twice the mean WIMP 
velocity) implies maximum momentum transfers of gmax — 550 MeV for nuclei with atomic 
weight A ~ 127. This q value is not small compared to the inverse nuclear size and one 
has to use the finite momentum transfer approximation for heavier nuclei. 

Despite the above-mentioned simple kinematic estimation of gmax (which is used over the 
text for illustration) it is worth noting that this gmax value is in fact too large and is almost 
not reachable. The WIMP-nucleus interaction is very weak, it occurs very rarely and therefore 
the impulse approximation for rather large g is well motivated and is used almost in all nuclear 
calculations reviewed in the paper. In the impulse approximation the WIMP-nuclear interac- 
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tion is described by means of the WIMP interaction with a nucleon from the nucleus A (see 
Appendix) and the maximal momentum transfer in the WIMP-nucleon system is considerably 
smaller than the gmax- Therefore for heavy enough nuclei, in general, the transfer momentum 

q < gmax- 

The full momentum dependence of the form factors must be calculated from detailed nuclear 
models, and the results are especially important for heavier nuclei 27|. Unfortunately, the 
sirnple phenomenological analysis in the OGM (odd group model) and EOGM (extended OGM) 
of [53] cannot be extended to the finite momentum transfer case, because the experimental data 
directly related to the neutralino-nucleus elastic scattering is not available Quite a number 
of high multipoles can now contribute, some of them getting contributions from components of 
the wave function which do not contribute in the static limit (i.e. at g = 0). Thus, in general, 
sophisticated Shell Model calculations are needed to account both for the observed retardation 
of the static spin matrix element and its correct dependence on transfer momentum, q. For 



the experimentally interesting nucl"'="^ oTrcfoTnc agq^ ^.n^ 73r 
been performed by Ressell et al. 



ear systems ^|Si and 32Ge very elaborate calculations have 
4ol |. In the case of 32 Ge a further improved calculation by 
Dimitrov, Engel and Pittel was carried out j3] by suitably mixing variationally determined 
triaxial Slatter determinants. Indeed, for this complex nucleus many multipoles contribute 
and the needed calculations involve techniques which are extremely sophisticated [4S]. Now 
the necessity for more detailed calculations especially for the spin-dependent component of the 
cross sections for heavy nuclei is well motivated. 

Further available sets of spin structure functions S^{q) are collected for nuclei from ^^F up 
to the 207pb. 

A. Fluorine, ^^F 

Fluorine- 19 is the isotope most sensitive to the spin-dependent WIMP- nuclear interaction 
j^l and a lot of experimental grou ps hope to explore this feature of ^^F experimentally (see. 



for example. 







15|, 



55 



56 



The finite momentum transfer spin structure function S^^{q) for ^^F, one of the lightest DM 
target medium, (together with other light targets ^^Na and ^^Si) was obtained for the first 
time by Vergados with co-authors [28]. The spin contribution to the differential cross section 
was carefully estimated by the shell-model calculations in the sd shell using the Wildenthal 
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interaction (see, for example, [58|,|59|), which was developed and tested over many years. This 
interaction is known to reproduce accurately many nuclear observables for sd shell nuclei. The 
Wildenthal two-body matrix elements as well as the single-particle energies are determined by 
fits to experimental data in nuclei from A = 17 to A = 39. The shell- model wave functions 
used by the authors were tested in the calculation of the low-energy spectra and ground state 
magnetic moment. Rather good agreement between the theoretical and experimental results 
was achieved (see Tables 2-4 in our previous revie w lZiI ). Their spin matrix elements are in 
good agreement with those of previous calculations j40|. 

The pure isoscalar, Sqq, pure isovector, S\\, and interference, SH terms (jH) of the fluorine 
(J = 1/2) structure function S^^{q) can be given in the form j28|: 

Sll{q) = X (2.610) X e- {pj,i) W + ^(li) ' 

2 7 _|_ 1 

Sll{q) = X (2.807) X {g?o,i)(«) + ^?(2,i)(«)} ' (6) 

SoUq) = (2.707) xe-«{P(o,i)(w)g(o,i)(«) + ^(2,i)(w)g(2,i)(w)}. 

Here the following functions are introduced: 

P(o^i)(m) = 0.1145^2 - 0.6667m + 1, g(o,i)(^^) = 0.1088^^ - 0.6667m + 1, 

P(2,i) (u) = -0.0026 + 0.0100 u, Q(2,i) (u) = 0.0006 u'^ + 0.0041 u. 

The dimensionless variable u = (g6)^/2 and the energy transfer to the nucleus A are con- 
nected by the relation E = ub'^/MA (or u = h'^EMA). Here h = l.OOA^/^ = 1.63 fm = 
1.63/0.1973 GeV^^ is the oscillator size parameter for ^^F. For Vj^ax ~ 600 km/s one has 
Mmax ~ 0.17 (or gmax ~ 71 McV, and -Emax ~ 140 keV) and Mmin ~ 0.011, which corresponds to 
the energy threshold of 10 keV (and q^^^ ^ 17 MeV). 

Following Ql we note that at a relatively large momentum transfer the nucleonic axial 

(c? ■ fyjp 

current gets modification (PCAC) a ^ a -. Therefore structure functions ® become 

more complicated: 

Smil) = ^ {Pfo.i)(«)(l + M^)) + ^(2,i)M(l + /32iu)) - 2/3o2(u)P(o,i)(«)P(2,i)(«)} , 
9 807 

Suii) = e-" {g?o,i)(«)(i + Po{u)) + g?2,i)(«)(i + P2{u)) - 2/?o2(«)g(o,i)(«)g(2,i)(«)} , 

SoUq) = ^e-"{P(o,i)(M)g(o,i)(w)(l + /?o(w)) + P(2,i)(w)g{2,i)(w)(l + /92(«)) (7) 

-/?02(m) X {P(^o,l){u)Q(2A){u) + Q{0,l){u)P(^2,l){u))} , 
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where 

Hu) = 
When m„ 

0.22 
0.2 
0.18 
I 0.16 



U + Ut^ 



u 



oo formulas, ^ pass into The fluorine structure functions calculated in 
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FIG. 1: 

(bottom), 

(top) for ""t as a 
the recoil energy E = 
in keV calculated by (0). With 
t'max ~ 600(700) km/s for the ^^F 
target one has E'max ~ 140(190) keV. 



(middle), 
^^F as a function of 

= ub-^/MA 



accordance with ((7j) are given in Fig. ^ as functions of the recoil energy. Further details can be 
found in the original paper 28 1. 



B. Sodium, "'•^Na 



In the modern most promising scintillator dark matter detectors (like, for example, the 



DAMA one) iodine is always used together wit 



(Nal) crystals (see for example, j3l| 



6G 



61 



6^ 



1 sodium in the form of large sodium iodide 
3)- The nucleus ^^Na (J = 3/2) lies in the 



41| can be 



middle of the sd shell and therefore the methods applied to other sd-shell nuclei 

n 

used. Ressell and Dean 4l|| have performed calculation for ^^Na exactly analogous to those 
done for ^^Si in |40| and for ^^Al in j42], including the use of harmonic oscillator wave functions. 
The following zero-transfer-limit spin structure parameters are obtained: (S^^) = 0.2477 and 
(S^^) = 0.0198 4^. The fits to the structure functions Sf^{q) as third order polynomials 



in y = {qb/2Y are given as follows 4l|: 



sll{y) 
sll{y) 
Sll{y) 



0.0380 - 0.1743 y + 0.3783 - 0.3430 
0.0647 - 0.3503 y + 0.9100 - 0.9858 ?/^ 
0.0275 - 0.1696?/ + 0.5077 - 0.6180 y^. 



(9) 
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Here y is obviously well below 1, the functions are rather accurate to y < ymax = 0.1875 
{Qmax ~ 100 MeV), which corresponds to the maximum halo velocity of f max = 700 km/ s and the 
heavy enough WIMP mass. Here the oscillator parameter b = 1.6864 fm = (1/117.01) MeV~^ 
is used. 

Using their shell-model ajjproach Vergados with co-authors also calculated the spin structure 
function S'^^{q) for '^^Na ^]. Their three terms of the S'^^{q) (J = 3/2) are 

2 7 + 1 

S^M = X (0.478) X e-" {P(li)(«) + P(li)(«) + ^(l,3)(«) + ^(l3)(«)} , (10) 

S^Kq) = X (0.346) X {g?o,i)(«) + ^?{2,i)(«) + Q{2,3)(«) + Q(4,3)(«)} , 

2, J + 1 

SoUq) = (0.406) xe-"{P(o,i)(n)g(o,i)(M) + P(2,i)(^/)Q(2,i)(u) + 

+P{2,3)(m)Q{2,3)(^) + P{A,3){u)Q{4,3){u)} , 

where P(o,i) (u) = 0.0477w^ - 0.6667m + 1, Q(o,i) (u) = 0.0465m^ - 0.6667m + 1, 

P(2,i)(m) = -0.0177 M^ + 0.1048m, Q(2,i)(m) = -0.0349^^ + 0.1494m, 

P(2,3) (m) = -0.0767 M^ + 0.6092 m, Q(2,3) (m) = -0.0894 + 0.7405 m, 

P(4 3) (m) = 0.0221 M^, g(4^3) (m) = 0.0287 mI 

The dimensionless variable u = {qhY/2 = 2y, b = l.OOA^^^ = 1.69 fm is the oscillator size 
parameter for ^^Na. Relevant spin static characteristics of j3l (S^^) = 0.2477 and (S^^) = 
0.0199 are almost exactly equal to ones from ,41]. In Fig. |2lthe recoil energy dependence of 
the spin structure functions calculated for ^^Na by © and (fTUjl is presented. One can see from 
Fig. El that the isoscalar spin functions S^q have practically the same behavior in the left and 
right panels, but, despite the same normalization at g = 0, the isovector Sff and mixed SqI 



spin functions from Vergados et al. [28] (in the right panel) are systematically a bit larger than 



ntpc 

m I 



the Sff and S^f functions from Ressell and Dean |41[ (in the left panel of Fig. 

C. Aluminum, ^^Al 

The ^^Al nucleus is one of the active ingredients in a very high-resolution and low-threshold 



64| . Engel, Ressell, Towner, and Ormand in j4 



sapphire-crystal (AI2O3) CRESST detector 

carried out calculations with the Lanczos m-scheme shell-model code CRUNCHER |65!I. The 
m-scheme basis for ^^Al contains 80115 Slater determinants. The agreement of the calculated 
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FIG. 2: Spin structure functions for ^^Na 5*00 (middle), 5?? (bottom), and 5gf (top) versus the recoil 



energy. Left: Sf f ap proximations from Ressell and Dean |4ll | in accordance with Q. Right: Sff from 



Vergados et al. [2^ following With 



700 km/s for Na target one has E„ 



230 keV. 



spectrum with that measured in Al is very good. A similar calculation for Si was fulfilled 
by Ressell et al. 4^. In the shell model, the expectation value of any one-body operator 
and therefore the structure functions can be easily calculated. To evaluate the S^'^{q) Engel, 
Ressell, Towner, and Ormand used a nuclear mean-field potential which is closer to the Woods- 
Saxon rather than the harmonic oscillator form. The length parameter b = 1.73 fm for the 



oscillator functions and the standard parameters for the Woods-Saxon potential were used 42 
in calculation of the structure functions. The most accurate ^^Al structure functions of 42 
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FIG. 3: Structure functions 5go (mid- 
dle), Sfl (bottom), and S^l (top) for 
^^Al as a function of the recoil en- 
ergy calculated by (|TT|) . For ^'^Al 
and the WIMP maximal velocity of 
700 km/s the maximum momentum 
transfer is Qmax ~ 

117 MeV and 

^max « 270 keV. 



11 



with the Woods-Saxon single-particle wave functions can be given to a high accuracy by the 
third-order polynomials 

S^Kq) = 0.0930 - 0.4721?/+ 1.0600 1.0115 ?/^ 

S^Kq) = 0.0657 - 0.4498 y+ 1.3504 - 1.6851 ?/^ (11) 

S^l{q) = 0.1563 -0.9360?/ + 2.4578 y^- 2.7262 

where y = (6g/2)^. Figure El shows these functions SfJ for aluminum versus the recoil energy 
(see formulas (fTTjl ). These three functions allow one to calculate the spin-dependent cross 



42|. 



sections for the neutralino interaction with the ^^Al target 

D. Silicon, "^^Si 

The first accurate calculation of the g-dependence of the spin structure functions for WIMP 
scattering off ^^Si was carried out by Ressell et al. [4^. With a reasonable two-body interaction 
Hamiltonian Ressell et al. calculated nuclear wave functions in an appropriate model space. 
The accuracy of those wave functions was checked by comparison of the calculated excited state 
energy spectrum, magnetic moments, and spectroscopic factors with the experimental observ- 
ables. The checked ground state wave functions were used for further calculations of static and 
more complicated finite momentum nuclear matrix elements involved in DM search calculations. 
Finite momentum transfer matrix elements and cross sections for the spin-dependent elastic 
scattering of neutralinos from ^^Si and ^^Ge 14^1 were evaluated in that shell-model scheme. 
Both these isotopes have a great number of configuration mixing and require very large model 
spaces. 



In particular, for evaluation of the wave 
interaction of Wildenthal (see, for example, jsS 



"unctions for Si the universal sd shell-model 



,|59|) in a full sd shell- model space was used 







and calculations were carried out with the Lanczos m-scheme shell-model code CRUNCHER 
As in the case of ^"^Al, the m-scheme basis for ^^Si contains 80115 Slater determinants. 
For static spin matrix elements the following values were obtained: (S^^) = —0.002 and (S^^) = 
0.13 (given in Table 4 of Q). 

Figure |^ (left panel) shows the recoil energy dependence of the spin partial structure func- 
tions 5*00, Sff and Sgi for ^^Si originally calculated in [4^ as a function of y = (6g/2)^ with the 
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oscillator parameter b = 1.75 fm for ^^Si. With fmax = 600(700) km/s the maximal momentum 
transfer is g^ax = 0.108(126) GeV, y^ax = 0.23(0.31) and E^^^ = 213(290) keV. 

For user's convenience, in their earlier work Ressell et al. obtained rather simple pa- 
rameterizations of the full spin structure function S'^^(g) (quenching is included) 



Slliy) = 0.00818 (a^e-^-^^s, ^ iMaje-^-''^^ - 2.06aoaie-^-^^=^^) 



(12) 



which provides a highly accurate fit to S'^^{y) for y < 0.15. For y > 0.15 this parameterization 
begins seriously underestimate the true value of 5*^^(7/) [40j. The lack of validity of (fT^ at high 
y (or q) is not a serious problem for WIMP masses up to about 600 GeV/c^. Anyway, when 
ao ~ fli, parameterization (fT^ does not reproduce accurately the full result of |4i| for all y 
given in the left panel of Fig. |3 
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FIG. 4: Spin structure functions Sqq (top) Slf (middle), and Sqi (bottom) for ^^Si as a function of 
the recoil energy. Left: results of Ressell et al. |4fl]. Right: these structure functions from Vergados 
et al. Q following equations (UD). With v^^x. ~ 700 km/s, for ^^Si target one has E'max ~ 290 keV. 



Vergados with co-authors also calculated the spin structure function S'^^(g) for ^^Si in their 
approach of [28]. The pure isoscalar, Sqq, pure isovector, Sff, and interference, S^f terms of 
the silicon structure function Q S^^{q) {J'" = 1/2"^) are given in the form 

2J 



SooiQ) = X (0.208) X e-« {P(li)(«) + P(li)(«)} , 

2 7-1-1 

Snil) = (0.220) xe-"{Q?o,y(«) + g?2,i)M} 



167r 
2J+ 1 

87r 



(13) 



X (-0.214) X e-"{P(o,i)(^/)Q(o,i)(u) + P(2,i)(M)Q(2,i)(n)} 



13 



where P(o,i)(m) = 0.2843 - 0.6667m + 1, <5(o,i)(w) = 0.2710 - 0.6667m + 1, 

P(2,i) (n) = -0.0567 + 0.4566 m, Q(2^) [u) = -0.0621 + 0.4680 u. 

Here u = (g6)^/2 = 2?/. For a rather large momentum transfer, in full analogy with ^^F (see, 
formulas ((Tj)), one has 

SooiQ) = ^ {Pli)iu)il + /3o(tx)) + P(li)(«)(l + /?2(tx)) - 2/?o2(tx)P(o,i)(tx)P(2,i)(tx)} , 

SlUq) = ^ e-" {Q?o,i)H(1 + /^oM) + Q(2,i)(«)(l + f32{u)) - 2(3o2{u) g(o,i)(«)Q(2,i)(n)} , 

SoUq) = ^^e-"{P(o,i)(n)Q(o,i)(«)(l + /3o(«)) + P(2,i)(n)g(2,i)(«)(l + /32(w)) 

-/5o2(«) (P(o,i)(M)Q(2,i)(n) + g(o,i)(^^)P(2,i)(^^))} , (14) 

where Po{u), P2{u), and /5o2(m) are defined in (jHl). Figure HI (right panel) presents structure 



functions for Si from Vergados et al. 



28| following formulas ()14|) . Both sets of spin structure 
functions for ^^Si given in the left j40| and right j3| panels of Fig. EJare very similar. 

E. Potassium, ^^K 

In the case of ^^K the shell-model diagonalization needed for the calculation of the nuclear 
spin matrix elements requires severe truncations to the active model space. The problem is that 
^^K is near the boundary between the sd and pf shells and excitations of particles into higher 
shells can have significant effects that are often not well simulated by effective operators. Thus, 



42 1 used an alternative scheme based on 



for this nucleus Engel, Ressell, Towner and Ormand 
perturbation theory for the evaluation of spin matrix elements. The authors considered two 
different residual interactions. One is related to the one-boson-exchange potential of the Bonn 
type, but it is limited only to four or five important meson exchanges. The resulting interaction 
has a weak tensor-force component typical of Bonn potentials. The other is represented by full 
G-matrix elements of the Paris potential parameterized in terms of sums over Yukawa functions 
of various ranges and strengths. This interaction exhibits a strong tensor force. The quality 
of the wave functions obtained was judged in terms of magnetic moments and Gamow- Teller 
matrix elements, including meson-exchange currents, isobar currents, and other relativistic 
effects. The magnetic moments calculated with the help of both interactions differ only slightly 
from each other and showed good agreement with the corresponding experimental values. The 
same nuclear wave functions of '^^K were also used for the calculation of (S^^) and (S^^) and 
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considered in 



the structure function S (q) [42]. Three different sets of spin structure functions for K were 



42j |. They are the "single-particle" functions, the full functions obtained from 



the modified Bonn interaction, and the full functions obtained from the Paris-based G-matrix. 
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FIG. 5: ^^K spin structure functions Sqq 
(top), Sff (middle), and S^f (bottom) 
versus the recoil energy, calculated by 
()15|1 from Engel, Ressell, Towner and 
Ormand 0. With i^max = 700 km/s, 
for one has g^ax = 169 MeV, 

^max = 390 keV. 



There is rather strong reduction of Sfj in comparison with their single-particle values. The 
strongest reduction is in S^Q^q), which is reduced to 25% and 20% of the single-particle value for 
the two residual interactions. The preferred choice (corresponding to the Paris-based G-matrix 
approach) can be rather accurately reproduced by the following fourth-order polynomials in 
y = (6g/2)2 (with b = 1.84 fm for ^^K): 



^39(g) = 0.0095 -0.0620?/ + 0.1628 7/2- 0. 1943?/=^ + 0.0891 ?/^ 
S^l{q) = 0.0298 -0.2176?/ + 0.6236 1/2- 0.8144 + 0.4050 1/^ 
SiKq) = 0.0332 - 0.2319?/ + 0.6385 ?/2_ 0.7985?/^ + 0.3810 ?/^ 



(15) 



The spin structure functions ()15|) for '^^K are given in Fig. El These structure functions allow 
one to interpret experiments that can look for tracks due to the interaction of dark-matter 
particles with nuclei in ancient mica 

F. Germanium, ^^Ge 

Germanium isotopes, especially large-spin (J = 9/2) ^'^Ge, are considered to be the most 
promising material for a real long-running experiment aimed at direct dark matter search. The 
first experiment with a pure (enriched) ''^Ge target was successfully performed at Gran Sasso 



40| (together with ^^Si). An equally 
3|. These authors 
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by the HDMS collaboration [19|,|66|. However, there are fundamental difficulties in describing 
the spin content of ''^Ge due to its complicated collective structure. 

The ffist accurate calculation of the g-dependence of the spin structure functions for WIMP 
scattering off '^'^Ge was carried out by Ressell at al. in 
comprehensive calculation was realized by Dimitrov, Engel and Pittel 
obtained significantly different (and improved) results in comparison with other studies and 
they argue that their results are more reliable than the previous ones. 

For the study of '^'^Ge Ressell et al. 0] chose the Petrovich-McManus-Madsen-Atkinson 
interaction jo?!, which is a reasonable approximation to a full G- matrix calculation. This in- 
teraction proved to be both adequate and tractable in shell model applications. Two different 
model spaces were considered. The "small" space was determined by an m-scheme basis di- 
mension of 24731 Slater determinants. The "large" space allowed much more excitations with 
an m-scheme basis dimension of 117137 Slater determinants. Despite fairly large size of the 
bases, rather severe truncations in the space were implemented. The small space is the small- 
est one in which it is possible to obtain agreement with the experimental spectrum energy 
levels. The dimension of the large basis was limited by the computer time and the memory 
storing constraints . No phenomenological interaction has been developed for Ge-like nuclei 
and fairly severe truncations to the model space have to be imposed to obtain manageable 
dimensions. The large model space wave function of ^^Ge led to an improved description of 
the ground state expectation values, in particular of the value of the magnetic moment, in 
comparison with previous estimates. The calculated magnetic moment /i from exceeds 
the experimental value, but the authors stressed that the same quenching of both /i and the 
Gammov- Teller spin matrix elements was almost universally required in shell model calcula- 
tions of all heavy nuclei. Assuming the isovector spin quenching factor to be 0.833, agreement 
with the measured fi is obtained. In principle, it is not obvious that quenching is really needed 
in neutralino-''^Ge scattering but if so, Ressell et al. believed that the correct answer might 
)e in the range between the quenched and unquenched values. It was found (see Table 6 of 
3]) that the zero-momentum-transfer spin-neutron matrix element (S^'^) of ^^Ge was a factor 
of 2 larger than the previous predictions. Thus, even if quenching is assumed, the calculated 
scattering rate is about twice as large as any of the estimates obtained before. Therefore Ressell 
et al. predicted a higher sensitivity for germanium dark matter detectors 
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Figure El (left panel) shows the recoil energy dependence of the partial structure functions 
S'oo, Slf and S^f for '^^Ge calculated by Ressell et al. in in terms of ?/ = (6g/2)^. For ''^Ge 
the oscillator parameter b = 2.04 fm, the maximal momentum transfer (fmax = 600 km/s) is 
9max = 0.271 GeV and -Emax = 537 keV. As in the case of ^^Si (see expression dH))) Ressell 
et al. gave rather simple parameterizations of the full spin structure function S'''^(g) (where 
quenching is included and which is valid for y < 0.2): 

Slfiy) = 0.20313 {1. 10246-' -^"^^^ + a^e-s-s^ss/ _ 2.099aoaie-^-^^^y). (16) 

The lack of validity of (jl6|) for ''^Ge at high y (or q) is not a problem for WIMP masses up to 
about 600 GeV/c^. When ao ~ ai, parameterization (fTI)|) does not reproduce accurately the 
full result for all y [4^ . 

A different sophisticated approach to evaluation of the spin structure of ^^Ge was considered 
by Dimitrov, Engel and Pittel in 45|. It relies on the idea of mixing variationally determined 
Slater determinants, in which symmetries are broken but restored either before or after vari- 

n n 

ation. This approach is described in detail in |6^. In the calculation of |45| the symmetries 
broken in the intrinsic states are those associated with rotational invariance, parity, and axial 
shape. The hybrid procedure used restores axial symmetry, parity invariance, and approximate 
rotational invariance prior to the variation of each intrinsic state. Subsequently, before mix- 
ing the intrinsic states the rotational invariance is fully restored. The procedure allows f ully 
triaxial Slater determinants at the expense of particle-number breaking. The results of 63| 
indicate that the trading of number nonconservation for triaxiality is a good idea, despite the 
apparent loss of pairing correlations traditionally associated with the former. Pairing forces 
evidently induce effective triaxiality. The numerical results j6^ show that the approach is ac- 
curate and efficient for describing even-even systems while also providing reliable reproduction 
of the collective dynamics of odd-mass systems 

For ''^Ge the calculations in 4^ were performed by assuming, both for protons and neutrons, 
a single-particle model consisting of the full Of,lp shell and the OgQ/2 and 0(^7/2 levels. The 
main idea was to include all of the single-particle orbits that could play an important role in 
reproducing low-energy properties of ^^Ge It is well-known that a crucial ingredient in 

any realistic nuclear-structure calculation is the appropriate form of the nuclear Hamiltonian. 
The one- and two-body parts of the Hamiltonian have to be compatible with each other as 
well as with the model space. This is difficult to achieve because microscopic two-body in- 
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FIG. 6: Spin structure functions Sqq (top) S'Jf (middle) , and S^f (bottom) for '^^Ge as a function of 
the recoil energy. Left: results of Ressell et al. l40ll . Right: the same structure functions, but from the 
"hybrid" method of Dimitrov, Engel and Pittel |45l | following equations ()17|) . With Vmax 

PS 600 km/s, 

for the '^^Ge target one has E'max ~ 537 keV. 

teractions, derived for example from a G-matrix, include monopole pieces that are unable to 
describe the movement of spherical single-particle levels as one passes from the beginning to 
the end of a shell. A procedure proposed for avoiding this problem consists basically in remov- 
ing all monopole components from the two-body interaction and shifting their effects to the 
single-particle energies. This procedure was used by Dimitrov, Engel and Pittel 45| — their 
two-body force was a fit to the Paris-potential G-matrix modified as described above. The 
calculated ground-state magnetic dipole moment is in good agreement with the experimental 
value. Ressell et al. in their large-space shell-model calculation were able to reduce fi 
significantly to — 1.239/iAr (with experimental value — 0.879/iAr) but could not account for the 
remaining difference. On the contrary, the calculation of Dimitrov, Engel and Pittel, despite 
the small number of intrinsic states, contains the full quenching required by the experimental 
data 45] . By making a comparison with the results of Ressell et al. 40] , significant disagree- 
ment is found for the neutron spin again. The calculated value of Dimitrov, Engel and Pittel 
is significantly smaller (— 0.920;Ujv)- The differences in the spin contributions, unlike those in 
the orbital angular momenta, strongly affect the WIMP-Ge scattering cross sections. Thus, 
contrary to Ressell et al. [4^, no significant increase can be expected in the neutralino-^'^Ge 
scattering rate in accordance with j45|. 
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45[ . Figure IHl (right panel) 
that determine the spin- 



The advantage of Dimitrov, Engel and Pittel's approach to calculation of neutralino cross 
sections is that it correctly represents the spin structure, requires neither quenching at g = 
nor arbitrary assumptions about the form factor behavior at g 7^ 
shows partial structure functions of Dimitrov, Engel and Pittel 4| 
dependent cross sections of elastic neutralino scattering off ^^Ge. 

Comparing the results for S^Q^q) (the pure isoscalar form factor) and Slf{q) (the isovector 
form factor) with the corresponding large-space results of Ressell et al. 4^ given in Fig. IHl (left 
panel) one can conclude that both 6*00 (V) and Slf{q) of Dimitrov, Engel and Pittel are reduced 
relative to the curves of Ressell et al. 40| . 

The polynomial fits, in terms ofy = (6g/2)^ with b = 2.04 fm being the oscillator parameter 
for ^^Ge, which well represent the structure functions in Fig. IHl (right panel), have the forms of 
the following sixth-order polynomials 

^73,^^) ^ 0.1606 - 1.1052?/ + 3.2320?/2 - 4.9245?/^ + 4.122%^ - 1.8016?/^ + 0.3211?/^ (17) 
Slf{y) = 0.1164 - 0.9228 y + 2.9753?/2 - 4.8709 y^ + 4.3099?/^ - 1.9661?/^ + 0.3624/, 
5o^3(y) = -0.2736 + 2.0374?/ - 6.2803 y^ + 9.9426?/^ - 8.5710?/^ + 3.8310?/^ - 0.6948 y^. 

Here we have to stress that expressions (fT7j) are valid only for ?/ < 1 (the smaller y the better 
the accuracy). At the same time for '^'^Ge, when the maximal WIMP velocity fmax = 600 km/s, 
one already has the maximal possible momentum transfer gmax = 0.271 GeV (-Emax = 537 keV, 
Qmaxb = 2.08) and yma.x = 1-96 > 1. Therefore formulas (fTTj) are obviously not valid over the 
full range of relevant momenta transfer. From Fig. IHl (right panel) it seems, perhaps, rather 
safe in practice to put all Sj^{y) = for ?/ > 1. 

Finally, it will be a very hard task to substantially improve the calculations of Ressell et 
al. [4^ and Dimitrov, Engel and Pittel j^l for these spin matrix elements. Nevertheless, 
Dimitrov, Engel and Pittel have proposed 45] several possible improvements to their analysis 
in the future. The most important is an explicit incorporation of the remaining orbits from 
the 2s,ld,0g shell, whose effects were treated very roughly in [45]. This should enable one to 
remove some of the arbitrariness in the single-particle energies that resulted from incomplete 
treatment of parity-mixing effects, etc (for details see original paper 0). 
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G. Niobium, ^^Nb 



The niobium isotope ^^Nb is an odd-group proton nucleus with a large ground-state angular 
momentum J = 9/2. It is a heavy enough nucleus, which can be represented by a basic shell- 
mode, space eo.xespo„di„g to three protons in the Ip.,. or 0,,, levels and two neutrons in 
the 1^5/2 level. This model space was considered by Engel, Pittel, Ormand and Vogel |46|]. In 
order to obtain better agreement with the experiment the authors extended it including in the 
"large" basis all states in which one proton or one neutron is excited from the above-mentioned 
"small" space to any level in the srf^f-shell. The resul ting space contains about 2700 states. 
They ended up with the nuclear spin matrix elements |46j (S^^) = 0.46 and (S^^) = 0.08 
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FIG. 7: Partial structure functions 
Slliq) (middle), Sll{q) (top) and 
Sll{q) (bottom) in ^SNb func- 
tion of the recoil energy obtained 
from paper of Engel, Pittel, Or- 



mand and Vogel 



Note that for 



^^Nb, when the maximal WIMP ve- 
locity Wmax = 600 km/s, one has 



0.345 GeV/c (and q, 



2 

max 



0.12 GeVVc^), En 
and ymax = 3.48. 



684 keV 



With relation (0) these values can be used to check normalization of the partial spin structure 
functions Sff{q) which are available only graphically from Fig. 3 of (and are given here in 
Table mi). In Fig. [7|we present recoil energy dependence of the Sff structure functions from 
Table US 



TABLE II: Tabulated partial spin structure functions Sff{q) from |4^ with in GeV'^/c 
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An interesting observation concerning interplay between spin-dependent and spin- 
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independent g-dependencies of nuclear structure functions is given by Engel in |46[. He noted 
that only nucleons near the Fermi surface contribute significantly to spin-dependent scattering. 
Because their orbits extend further out than those in the core, the form factor near g = 0, 
which reflects the mean square radius of the contributing nucleons, will always decrease faster 
for spin-dependent scattering than for spin-independent scattering. By the same token though, 
the Fermi-surface nucleons have a higher momentum on the average, and the spin-dependent 
form factors will therefore be larger at high q than their spin-independent con terp artners. 

On the other hand, for rather light isotopes (F, Na, Si) Vergados et al. 2S| claimed the 
opposite — the drop of the scalar form factor with increasing q is less dramatic compared to 
that of the spin structure functions. Therefore, contrary to heavy-mass targets, the light-mass 
targets have a better sensitivity to the scalar interaction for rather large q than to the spin 
one. Nevertheless, this advantage looks illusive due to the weak A^-enhancement of scalar 
interactions for the low-mass isotopes. 

H. Tellurium, I23,i25rpg 

The Theory of Finite Fermi Systems was used by Nikolaev and Klapdor-Kleingrothaus 
to define the g-dependence of nuclear form factors (structure functions) in spin-dependent 
WIMP scattering off the ^^^Te isotope. This approach allows one to describe the reduction 
of single-particle spin-dependent matrix elements in the nuclear medium. Unfortunately, only 
structure functions for nucleus spin interaction with the pure bino-like {a^/ai = 0.297) lightest 
neutralino are presented in 0]. Therefore it is not possible to extract the partial structure 
functions S}^^{q) for the ^^^Te isotopes from 3]. 

Ressell and Dean have performed most accurate nuclear calculations of the neutralino- 
nucleus spin-dependent cross section for ^^^Te (together with ^^^'^^^Xe and ^^"^I) in 0]. The 
details of the shell model and residual nucleon-nucleon interactions used in the calculations 
are given in the next Section. The g-dependence of complete spin structure functions S}^^{q) 



for ^^^Te one can br found in Fig. 2 of 4l|. Ressell and Dean gave rather simple parame- 

terizations of the complete structure functions Sl^^{q) as tables of the coefficients Ck (given 

6 

in Table Hn)) of 6th order polynomials in y: Slj^{q) = '^C^y'', where y = (g6/2)^ and 

A:=0 

b = 2.24 fm = 11.35/GeV for ^^^Te. These so-called abbreviated structure functions are only 
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valid for y < 1 (or g = < 176 MeV) and therefore for the recoil energy E < 124 keV. 



TABLE III: ^^^Te isotope. The first column gives the order of y'^, the next three columns give the 
corresponding values of the for Sqq^, -Sqi^' ^11^ for the Bonn A calculation. The last three 
columns present the same results for the Nijmegen II calculation. Prom 41 1. 
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Nijmegen II 




0125 
'-'00 


q125 
'-'01 


q125 
'-'11 
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1.9020 


-1.0944 
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2.3905 
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0.2380 
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-0.5245 


0.3241 



TABLE IV: ^^^Te isotope. The first column gives the order of y^, the next three columns give the 
corresponding values of the Ck for Sqq^, -Sqi^' "^iP fo'^ Bonn A calculation. The last three 
columns present the same results for the Nijmegen II calculation. From 41 1. 
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Ressell and Dean also gave their best parameterization of the complete structure functions 



Slf^{q) in the form of 8th order polynomials in y multiplied by a factor of exp{—2y): 



■2?; 



;i8) 



All coefficients Ck for the Te isotope are given in Table IIVI These so-called full structure 
functions S}J^{q) are valid for all kinetically available y < 10 (or q < 557 MeV) and are 
presented in Fig. |H1 for the Born A potential (three left columns in Table IIV|) . 



22 




1 mvJ^'Mimnk- 



■ 
I 



100 200 300 

Recoil Energy, keV 



400 



FIG. 8: Partial structure functions 
SooHq) (top), S^f{q) (bottom) and 
Slf^iq) (middle) in ^^^Te func- 
tion of recoil energy obtained from 
Ressell and Dean j4l| by means of 
parameterization (|18jl. Note that 
for ^^^Te, when the maximal WIMP 
velocity fmax = 600 km/s, one has 
gmax = 464 MeV/c, S^nax = 919 keV and 

ymax = 6.95. 



I. Iodine, ^'^'^I 

The iodine isotope ^^^I is a decisive component of large sodium iodide (Nal) crystals and 
therefore plays the most important role in the modern dark matter search (see for example 



DAMA |3l|,|60| and ELEGANT, NAIAD and ANAIS experiments |6l|, 162, |63|). 

Ressell and Dean have performed the most accurate nuclear shell-model calculations of the 
neutralino- nucleus spin-dependent cross section for ^^^I in j4l!|. Within the framework of their 
approach Ressell and Dean considered two residual nuclear interactions based upon recently 
developed realistic nucleon-nucleon Bonn A [3| and Nijmegen II jl^ potentials. These two 
nucleon-nucleon potentials were used in order to investigate the sensitivity of the results to 
the particular nuclear Hamiltonian. The Bonn-A-based Hamiltonian was derived for the model 
space consisting of the 1(77/2, 2^5/2, 3si/2, 2^3/2, and l/iii/2 orbitals, allowing one to include all 
relevant correlations. In order to get good agreement with observables for nuclei with A ^ 
130, the single-particle energies (SPEs) were adjusted. The SPEs were varied until reasonable 
agreement between calculation and experiment was found for the magnetic moment, the low- 
lying excited state energy spectrum, and the quadrupole moment of ^^^I. The similar procedure 
was used by Ressell et al. in [4^. Once the SPEs are specified, a reasonable Hamiltonian can 
be used for the nuclei under investigation (^^''I, ^^^'^^^Xe and ^^^Te). The same scheme was used 
for the Nijmegen Il-based Hamiltonian. 

To perform a full basis calculation of the ^^^I ground state properties in the space consisting 
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of the 1(77/2, 2(i5/2, 3si/2, 2^3/2, and 1/^11/2 orbitals, one would need basis states consisting of 
roughly 1.3 x 10^ Slater Determinants (SDs). Current calculations can diagonalize matrices with 
basis dimensions in the range 1-2x10'' SDs. Therefore clearly severe truncations of the model 



space are needed 
truncation scheme that includes t 
(after relevant truncations, see 



4l|. Fortunately, given the size of the model spaces that can be treated, a 
le majority of relevant configurations can be devised. Finally 



4l| for details) the m-scheme dimension of the ^^^I model 



space is about 3 million SDs. The calculated observables agree well with experiment. These 
interactions do not seem to prefer excitation of more than one extra neutron pair to the l/iii/2- 
Most configurations have six neutrons in that orbital, while eight are allowed. Hence, this 
model space is more than adequate. It is this truncation scheme that was also used for the two 
xenon isotopes considered {A = 129 and 131). 



In almost every instance, the results of Ressell and Dean j41| (Tables 8-10 of [7]) show 
that the spin (Sj) (z = p, n) carried by the unpaired nucleon is greater than that found in the 
other nuclear models. Despite these larger values for (Sj), these calculations have the magnetic 



moment in good agreement with experiment in all cases. The larger values o: 



(Sj) are due to the 

fact that more excitations of the even group of the nuclei were allowed j4l|. There are visible 
differences in the response due to the two different forces. In all cases reasonable agreement 
between calculation and experiment for the magnetic moment (using free particle ^^-factors) is 
achieved. It is obvious that the differences between the two calculations are non-trivial but 
they are quite a bit smaller than the differences coming from the use of alternative nuclear 
models. This shows that the interaction is not the primary uncertainty in calculations of the 



neutralino-nucleus spin cross sections 



4l|. The results obtained by Ressell and Dean give a 



factor of 20 increase in iodine's sensitivity to spin-dependent scattering over that previously 
assumed. Due to the form factor suppression the spin response of the sodium iodide detector's 
is still dominated by ^^Na but not to the extent previously thought 0|. 

The reduced matrix elements of the multipoles in the definition of the S^{q) (jH)) are easily 
evaluated in the harmonic oscillator basis in the nuclear shell model. Almost all calculations 
of S^{q) used bases of these harmonic oscillator wave functions. In Ressell and Dean used 
more realistic Woods-Saxon wave functions to evaluate 0. To check influence of the wave 
function basis the Bonn A structure function S™{q) for the ™\ isotope was also calculated 
with the harmonic oscillator wave functions. The Woods-Saxon wave functions made a signifi- 
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cant difference at extremely high momentum transfers when compared to the usual harmonic 
oscillator wave functions. At more modest momentum transfers the difference is found to be 
small. 

The oscillator parameter, b = IfmA^/^, is usually retained as the size parameter in 



Woods-Saxon evaluations of 5'^(g). In |4l| a standard S(i-shell parameterization is used: 
b = (41.467//iu;)^/2 ^j^^^ = 45^-^/3 _ 25v4~2/3 j^^y^ Therefore for all isotopes with 

A ^ 127 one has b = 2.282 fm. In general, for heavy enough nuclei (say, with A > 100) it is 
especially useful to present structure functions in terms of dimensionless variable y = (g6/2)^. 
For y <^ 1 {y > 1) the effects of finite momentum transfers are usually rather small (quite 
large). For these nuclei ?/max = (Q'max&/2)^ ~ 10 ^ 1, and nuclear structure form factors are 
very important. Nevertheless these extremely large values of y are only valid for extremely mas- 
sive WIMPs (m^ ^ m^) moving with an almost escape velocity (fmax ~ 700 km/s). A more 
realistic WIMP with mass of about 100 GeV/c^ moving at an average velocity f ~ (f ) = 10~^c 
would have t/max — 0.4. Anyway, in order to cover the entire elevant neutralino parameter 



space, the structure functions were evaluated up to y = 10 



4l|. 



The complete partial spin functions SfJy) for ^^''I, ^^^'^^^Xe and ^^^Te and y < 2 (g^ ^ 60000 



r2 



MeV ), originally calculated by Ressell and Dean with the Bonn- A- and Nijmegen-II-based 



Hamiltonians, are presented graphically in Fig. 2 of 4l|. For practical use Ressell and Dean 



gave the simple parameterization of ^^''I structure functions Slj'^{q) = ^ C^y'' as 6th order 

k=0 

polynomials in y with Ck coefficients from Table |3 These so-called abbreviated structure 
functions are only valid for y < 1. Ressell and Dean also gave so-called full structure functions, 
which are valid for y < 10 and reproduce the complete structure functions S^'^'^{q) as 8th order 
polynomials in y multiphed by a factor of exp(— 2?/): 

Sl^\q)=e-'yj:C,y>' (19) 

k=0 

All coefficients Ck for the ^^^I isotope are given in Table IVII For example, with the entries in 
the Bonn A section of Table IVll one obtains the analytical forms for ^^''I spin stricture functions: 



Sl^\y) = e-'^y (0.0983 - 0.4891|/ + 1.14022/^ _ 1.4717?/^ + 1.1717?/^ (20) 

-0.5646?/^ + 0.1583y^ - 0.0239?/^ + 0.0015?/^) , 
Slf{y) = (0. 1199 - 0.6184?/ + 1.5089?/2- 2.0737?/^ + 1.7731?/^ 
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-0.9036?/^ + 0.2600y*^ - 0.0387y^ + 0.002%^) , 
Slf{y) = (0.0366 -0.1950y + 0.5049?/2- 0.7475?/^ + 0.7043/ 
-0.3930?/^ + 0.1219/ - 0.0192/ + 0.0012/) . 



TABLE V: 



127 



I isotope. The first column gives the order of y'', the next three columns give the 



corresponding values of the Ck for Sqq'^ , S™ , and 5jp for the Bonn A calculation. The last three 



columns present the same results for the Nijmegen II calculation. From 41 1 





Bonn A 


Nijmegen II 




0127 
'-'00 


q127 
'-'01 


q127 


0127 
'-'DO 


q127 
'-'01 


q127 
'-'11 


/ 


0.0983 


0.1199 


0.0365 


0.1165 


0.1619 


0.0562 


/ 


-0.6750 


-0.8436 


-0.2627 


-0.7923 


-1.1403 


-0.4085 


/ 


2.1353 


2.7354 


0.8751 


2.4985 


3.7144 


1.3778 




-3.7595 


-4.9303 


-1.6146 


-4.3831 


-6.7158 


-2.5702 


/ 


3.7774 


5.0581 


1.6908 


4.3850 


6.8938 


2.7087 


yh 


-2.0091 


-2.7361 


-0.9302 


-2.3222 


-3.7259 


-1.4945 


y6 


0.4355 


0.6008 


0.2069 


0.5015 


0.8171 


0.3329 



™1 isotope. The first column gives the order of y^, the next three columns give the 
values of the for Sqq^, SqV, and for the Bonn A calculation. The last three 



TABLE VI: 
correspondin 

columns present the same results for the Nijmegen II calculation. Prom [41 1 





Bonn A 


Nijmegen II 


x(e-22/) 


ol27 
•^00 


Ql27 
'-'Ol 


ol27 
'-'ll 


Ql27 
'-'00 


ol27 
'-'Ol 


Ql27 
'-'ll 


/ 


0.0983 


0.1199 


0.0366 


0.1166 


0.1621 


0.0563 


/ 


-0.4891 


-0.6184 


-0.1950 


-0.5721 


-0.8363 


-0.3038 


/ 


1.1402 


1.5089 


0.5049 


1.3380 


2.0594 


0.7948 


yi 


-1.4717 


-2.0737 


-0.7475 


-1.7252 


-2.8319 


-1.1703 


/ 


1.1717 


1.7731 


0.7043 


1.3774 


2.3973 


1.0637 




-0.5646 


-0.9036 


-0.3930 


-0.6700 


-1.2121 


-0.5713 


y6 


0.1583 


0.2600 


0.1219 


0.1905 


0.3486 


0.1722 


/ 


-0.0239 


-0.0387 


-0.0192 


-0.0292 


-0.0522 


-0.0266 


y8 


0.0015 


0.0024 


0.0012 


0.0019 


0.0032 


0.0017 



The structure functions Sjj (y) for the I isotope and for both nucleon-nucleon forces 



(Bonn A and Nijmegen II) up to y = 10 were carefully analyzed in j4l|. Some similarities and 
differences between these two sets of iodine structure functions were observed. Despite that 
the differences between the two calculations are non-trivial but they are quite smaller than 
the differences coming from the use of alternative nuclear models. This shows, in particular, 
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FIG. 9: Partial structure functions 

of Ressell and Dean Sjlf'^ for ^^'^I 

''J 

as a function of recoil energy cal- 
culated from formulas ()2U() . Note, 
for ^^^I, when the maximal WIMP 
velocity fmax = 600 km/s, one has 
gmax « 472 MeV/c, ^^ax « 934 keV and 

ymax ~ 7.44. 



that the interaction is not the primary uncertainty in calculations of the neutralino-nucleus 
scattering cross sections Q|. 

Finally, it is perhaps a right place to note that for y ^ Umax ~ 7.5, which already corresponds 
to a rather large energy (-Emax ~ 900 keV) transmitted to the target nucleus, one may also can 
expect a non-negligible contribution from inelastic WIMP-nuclear interactions not considered 
here. 

J. Xenon, ^^ixe and ^'^^Xe 

Xenon (gaseous and liquid) is a very popular target material for modern large-scale dark 
matter detectors (see, for example, [t^jIi^Qi)- -^^^ ^^^^ time spin-dependent scattering of 
SUSY-like dark matter particles from ^^^Xe (J = 3/2) nuclei at non-zero momentum transfer 
was considered by Engel in Q . The configuration- mixing quasiparticle Tamm-DancofF approx- 
imation (QTDA) was used. In the zeroth order the ground state of ^^^Xe was represented as 
the 1(^3/2 quasineutron excitation of the even-even core |0) treated in the BCS approximation 
(BCS-based model of the Fermi surface). In order to incorporate nuclear structure correc- 
tions originating from the residual interaction three-quasiparticle configurations of the form 

['^d3/2['^l-'^/]^]^''^l'-') ^"^^ ['^d3/2Ki^/]'^]^^^l'-') were admixed. Here vr^ and represent the proton 
and neutron quasiparticle creation operators, K is an arbitrary intermediate angular momen- 
tum, and k, I run over a valence space consisting of the 2s, Id, Og and Oh harmonic oscillator 
levels jl I . Despite the fact that the amplitudes associated with the admixed three-quasiparticle 
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states are small (less than 5%), these admixtures can lead to a substantial effect. The ex- 
perimental value of the magnetic moment of ^^^Xe, which is about 0.69 was reproduced 
with an accuracy of 2%. The same approximation scheme results in (S^^^) = —0.041 and 
(S^^^) = —0.236. Figure Uni (left panel) shows Engel's partial structure functions S}^^ in ^^^Xe 
as a function of the recoil energy, recalculated from Table IVIII 



TABLE VII: Tabulated spin structure functions Sjf^{q) from Fig. 3 of f?]. 








0.0025 


0.005 


0.01 


0.015 


0.02 


0.025 


0.03 


0.04 


0.05 


0.06 




0.4 
0.020 
-0.056 


0.0215 
0.009 
-0.028 


0.014 
0.006 
-0.019 


0.01 

0.004 
-0.013 


0.009 
0.003 
-0.01 


0.008 
0.0027 
-0.009 


0.0075 
0.0025 
-0.008 


0.0066 
0.0023 
-0.007 


0.005 
0.0019 
-0.005 


0.0035 
0.0015 
-0.003 


0.0017 
0.001 
-0.001 



Engel noted again that the spin-dependent cross section in ^^^Xe falls as the momentum 
transfer increases, but more slowly than the spin-independent cross section. The spin-dependent 
efficiency is higher than the spin-independent one, being substantially higher for very heavy 
neutralinos. The relatively long tail of the spin-dependent structure functions is caused by 
nucleons near the Fermi surface, which do the bulk of the scattering. The core nucleons, which 
dominate the spin-independent cross section, contribute much less at large q. These are very 
general statements that should apply to other heavy nuclei as well jl|. 

The Theory of Finite Fermi Systems (TFFS) was also used by Nikolaev and Klapdor- 
Kleingrothaus j3] to define the g-dependence of nuclear structure functions (form factors) 
in spin-dependent WIMP scattering off the ^^^Xe isotope (as well as ^^'^Te). The quenching 
effect (at zero and relatively low q), due to reduction of single-particle spin-dependent matrix 
elements in the nuclear medium, and its disappearance at higher q was observed in The 
same g-behavior was also discussed before by Engel at al. l|, |46|. Nevertheless the shape of 
the S^^^ia) from Q| differs a bit (at intermediate q) from the one obtained with the oscillator 
)asis in Unfortunately, partial structure functions Slj^{q) for ^^^Xe are not available from 
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The most sophisticated shell-model treatment of the spin structure functions in ^^^Xe and 
^^^Xe (as well as in ^^''I and ^^^Te isotopes) was performed by Ressell and Dean in |41j. The 
details of the calculations are given in the previous Section. Their complete sets of the structure 
unctions S'^^^'^'^^(g) are valid for all relevant values of the momentum transfer (see Fig. 3 in 
4ll|). Ressell and Dean gave the simple parameterizations of the complete structure functions 
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FIG. 10: Partial structure functions SQQ^{q) (top), SQf^{q) (bottom) and Slf^{q) (middle) in ^^^Xe as 
a function of the recoil energy. Left: results of 1| from Table IVm Right: the full parameterizations of 
[4l| for the Bonn A potential (three left columns in Table HXj) . Note, that for ^^^Xe, when the maximal 
WIMP velocity Vmax = 600 km/s, one has (/max ~ 487 MeV/c, E^nax ~ 963 keV and ymax ~ 7.92. 



TABLE VIII: ^^g.isixe isotopes. The first column gives the order of y'^ , the next three columns give 
the corresponding values of the Ck for 5qo^'^^^, Sq^^'^"^^ , and Sl'^'^^^ for the Bonn A calculation. The 
last three columns present the same results for the Nijmegen II calculation. From j4l| . 





Bonn A 


Nijmegen II 




q129 
'-'00 


cl29 
'-'01 


0129 
'-'11 


Ql29 
'-'00 


q129 
'-'01 


q129 
'-'11 




0.0713 


-0.1216 


0.0518 


0.0465 


-0.0853 


0.0392 




-0.4804 


0.8745 


-0.3949 


-0.3138 


0.6150 


-0.2991 




1.4726 


-2.8317 


1.3433 


0.9656 


-1.9847 


1.0087 




-2.5323 


5.0922 


-2.5152 


-1.6666 


3.5496 


-1.8648 


/ 


2.4968 


-5.1976 


2.6480 


1.6477 


-3.6023 


1.9399 




-1.3071 


2.7924 


-1.4557 


-0.8642 


1.9257 


-1.0564 


yd 


0.2796 


-0.6088 


0.3228 


0.1851 


-0.4182 


0.2326 




clSl 
'-'00 


cl31 
'-'Ol 




clSl 
'-'00 


Qiai 

'-'Ol 




yO 


0.0296 


-0.0545 


0.0251 


0.0277 


-0.0497 


0.0223 


y' 


-0.1852 


0.3676 


-0.1812 


-0.1754 


0.3389 


-0.1627 


y' 


0.5934 


-1.1813 


0.5932 


0.5604 


-1.1002 


0.5427 


y3 


-1.0351 


2.0529 


-1.0389 


-0.9969 


1.9709 


-0.9892 


y' 


1.0049 


-1.9827 


1.0071 


1.0100 


-1.9996 


1.0150 


y5 


-0.5078 


0.9967 


-0.5071 


-0.5402 


1.0681 


-0.5459 


yd 


0.1037 


-0.2026 


0.1031 


0.1174 


-0.2316 


0.1189 
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^129,131 j^g,^ as 6th order polynomials in y with coefficients Ck in Table IVTm S^{q) = ^ CkU^. 

k=0 

These so-called abbreviated structure functions are only valid for y < 1. 

For the ^^^'^^^Xe isotopes (as previously for ^^'^I and ^^^Te) Ressell and Dean also presented 
their parameterization of the full structure functions in the following analytical form 

SfjiQ) = f E Cw' + C9T^) e-'y. (21) 
The form is also valid for y < 10. The relevant coefficients Ck are given in Table HXl Figure ITUl 

TABLE IX: I29,i3ixe isotope. The first column gives the order of y^, the next three columns give 
the corresponding values of the Ck for Sqq^'^^^ , Sq^^'^^^ , and S^^'^^^ for the Bonn A calculation. The 
last three columns present the same results for the Nijmegen II calculation. From j4]] |. 





Bonn A 


Nijmegen II 


X e " 


Cl29 
^00 


cl29 
^01 




cl29 
^00 


cl29 
^01 


'-^ii 




0.07132 


-0.12166 


-2.05825 


0.04649 


-0.08538 


-1.28214 




-0.34478 


0.64435 


1.80756 


-0.22551 


0.45343 


1.09276 




0.75590 


-1.52732 


-1.27746 


0.49905 


-1.06546 


-0.71295 


y3 


-0.93345 


2.02061 


0.65459 


-0.62244 


1.38670 


0.31489 


/ 


0.69006 


-1.57689 


-0.22197 


0.46361 


-1.05940 


-0.08351 


y5 


-0.30248 


0.72398 


0.04546 


-0.20375 


0.47576 


0.01059 


y' 


0.07653 


-0.19040 


-0.00427 


0.05109 


-0.12208 


0.00023 


y' 


-0.01032 


0.02638 


-0.00014 


-0.00671 


0.01643 


-0.00024 


y' 


0.00057 


-0.00149 


0.00004 


0.00036 


-0.00089 


0.00002 


1 

l+y 


0.0 


0.0 


2.11016 


0.0 


0.0 


1.32136 




cl31 


cl31 
'-'Ol 


0131 
Oil 


cl31 


cl31 
'-'Ol 


0131 
Oil 


yO 


0.02964 


-0.05455 


0.02510 


0.02773 


-0.04978 


0.02234 


y' 


-0.13343 


0.27176 


-0.13772 


-0.12449 


0.24725 


-0.12206 


y' 


0.37799 


-0.72302 


0.36661 


0.32829 


-0.63231 


0.31949 


y' 


-0.57961 


1.05450 


-0.53851 


-0.48140 


0.89642 


-0.46695 


y' 


0.57890 


-0.97133 


0.49255 


0.47565 


-0.81645 


0.42877 


y' 


-0.34556 


0.53842 


-0.26990 


-0.28518 


0.45235 


-0.23679 


y' 


0.11595 


-0.16899 


0.08369 


0.09682 


-0.14267 


0.07408 


y' 


-0.02012 


0.02742 


-0.01340 


-0.01710 


0.02335 


-0.01197 


y8 


0.00142 


-0.00181 


0.00087 


0.00124 


-0.00156 


0.00079 



(right panel) presents parameterizations S}^^ from ()21|1 for the Bonn A potential (three left 
columns in Table HX]) as a function of the recoil energy. From Fig. El where parameterizations 
of the full structure functions Sff^ are depicted (see Table HXl for ^^^Xe), one can make a general 
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FIG. 11: Partial structure functions S^^^{q) (top), S^f{q) (bottom) and (g) (middle) in ^^^Xe 



from 4l| as a function of the recoil energy. Left: the full parameterizations for the Bonn A potential 
(three left columns in Table llX|) . Right: the full parameterizations for the Nijmegen II potential (three 
right columns in Table ITXj) . 



conclusion that there is no large difference between the calculations performed in j41| with the 
Bonn-A-based and Nijmegen-II-based Hamiltonians. 

Structure functions S^^^{q) for ^^^Xe have been calculated in the context of two other nuclear 
models mentioned above, the QTDA by Engel jl| and the TFFS by Nikolaev and Klapdor- 
Kleingrothaus [g^. Following Ressell and Dean ^41], we briefly touch upon the differences 
in S^^^{q) that are the result of using different nuclear models. All three calculations show 
significant quenching compared to the single-particle estimate. The spin distributions obtained 
in the QTDA and TFFS are somewhat different (see Table 10 from 7])) while the full structure 
functions (for bino neutralinos) are quite similar. While the values for (S„) differ very little 
between Ressell and Dean's and Engel's result in the QTDA, the difference in the values of 
5'^'^^(0) is almost a factor of 2 between the two calculations. It should be noted as well that 
both the QTDA and the TFFS calculations of S^^^{q) asymptotically reach the single-particle 
structure function. This is not the case in Ressell and Dean's calculations, which are well 
below the single-particle estimate for all values of q^. It is apparent that the shell-model- 



4l| . This difference 



derived structure functions have a much steeper fall-off as a function of 
can be considered as an estimate of the uncertainty followed from the choice of the nuclear 
model. 
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Finally, we note that Ressell and Dean in 41| gave accurate structure function parameteriza- 
tions for the ^^^'^'^^Xe, ^^''I and ^^^Te isotopes which are very useful in rather precise calculations 
of the event rates in heavy-target dark matter detectors. The calculations of Ressell and Dean 
contain more excitations within the model space and use more modern and realistic nuclear 
interactions than others in the literature, with a possible exception of the QTDA calculations 
Q for i3iXe lij. 



K. Lead, ^^'^Ph 



Among the nuclei which can be considered as targets for direct dark matter detection, ^°^Pb 
seems a potentially rather interesting candidate. The spin matrix element of this nucleus has 
not been evaluated quite accurately since one expected that WIMP-nuclear spin interaction is 
important only with light nuclei. But the spin matrix element in the light systems is quenched. 
On the other hand, the spin matrix element of ^''^Pb (J = 1/2), especially the isoscalar one, 
does not suffer unusually large quenching, as is known from the study of the magnetic moment 
It is believed that ^°^Pb has a quite simple structure, its ground state can be described 
as a 2pi/2 neutron hole outside the doubly magic (closed-shell) nucleus ^°^Pb. Due to its 
low angular momentum, only two multipoles L = and L = 2 can contribute even at large 
momentum transfers. This is why Kosmas and Vergados in 4^ chose ^°^Pb for investigation of 
momentum transfer dependence of spin matrix elements in the heaviest possible nuclei relevant 
to dark matter search. In the g = limit Vergados and Kosmas gave the spin matrix element 
in the simple form | JP = lf]lfio(0) + f\ni{0)\\ and found that l]o(0) = -0.95659/v^ and 
f2i(0) = 0.83296/v^ |4J, |4S|. The momentum transfer dependence of the total cross section 
for elastic scattering of cold dark matter candidates, i.e. the lightest su persy mmetric particle 
(LSP), from ^'''^Pb was examined by Kosmas and Vergados in j48 
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If the spin-dependent differential cross section is taken in the form of Engel, Ressell et al. 



4l| and the spin structure function h{q) has the form of Vergados et al. j28|: 

1 r. 



d(7 



Qfi2 

-h{q)dq\ Kq) 



2J+1 ' ' ^-4 
then the partial structure functions for ^°''Pb S^Q'^^q), Sff^{q), and S^f^^q) have the form |2^ 

2J + 1 



(0.305) /oo(g). 
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2J+ 1 

IGvr 
2J+ 1 

Stt 



(0.231) 

(-0.266) /oi(g). 



(22) 



The spin-dependent neutralino-nucleon couplings /^'^ are analogous to the parameters ao.i from 
(Q. Finite momentum dependence of these structure functions is concentrated in Iij{q), which 
can be defined, following |28l. l48 1 . as an integral over the forward scattering angle ^ = Pi • q ^ 



(23) 



with a(q) = (2J+l)-^/'(J||E^(j>«(j>"'''"^V) and uJo{j) = l, uj,{j) = r,{j). 

Here <j{j), T^i^j), are the spin, the third component of the isospin (r^lp) = \p)), and the 
coordinate of the jth nucleon, q is the momentum transferred to the nucleus. 

To a good approximation, the ground state of the ^82Pb nucleus can be described as a 2pi/2 
neutron hole in the ^82Pb closed shell ^|. Then for the L = multipole one finds 



Oi(q) = (l/V3)F2p(q') = -iloiq) and Jqo = loi = hi = 2 [\ [F2piq^)f = [F2piq^)f 

•J 

Here the form factor of a single-particle harmonic oscillator wave function has the form 



iVm ax 



u) 



(24) 



A=0 



with the dimensionless variable u = {b = 2.434 fm for Pb) and the coefficients 7^"'^ given 

in Table IXl With entries [n, I = 2,1) of Table 1X1 the form factor squared 



5 17 . 31 
1 u-\ 



9 4 1 
■ u — 



U'' 



' " V 3 " ' 15 " 105 ' 280 840 
defines the partial spin structure functions for ^°^Pb in the simplest case (pure 2pi/2 neutron 
hole in the closed shell and L = 0). 

Even though the probability of finding a pure 2pi/2 neutron hole in the | ground state of 
^°''Pb is greater than 95%, the ground state magnetic moment is quenched due to the 1"*" p-h 
excitation involving the spin orbit partners. Hence, Kosmas and Vergados expected a similar 



suppression of the isovector spin matrix elements 
has 



48|. In this case for the L = multipole one 



^j=o,i(q) 



Cl {F2p(g2)/v^ - 8 [(7/13)i/2CiFo.(g2) + (-l)^' (5/11)^/2^2^0,(9^)] }. 
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TABLE X: The coefficients 'Jx ^ entering into the polynomial describing the form factor l\24\i of a 
single particle harmonic oscillator wave function up to 6huj. From 
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-1/6 
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-1/3 


1/24 
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-1/3 
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1 1 


1 


-1/2 
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-1/2 


1/20 


-1/840 
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-2/3 


11/60 


-1/60 


1/1920 






1 2 


1 


-2/3 


19/120 


-11/840 


1/3360 






4 


1 


-2/3 


1/10 


-1/210 


1/15120 






2 1 


1 


-5/6 


17/60 


-31/840 


9/4480 


-1/26880 




1 3 


1 


-5/6 


29/120 


-47/1680 


37/30240 


-1/60480 
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-5/6 


1/6 


-1/84 


1/3024 


-1/332640 




3 


1 


-1 


17/40 


-31/420 


27/4480 


-1/4480 


1/322560 


2 2 


1 


-1 


2/5 


-1/15 


41/8064 


-1/5760 


1/483840 


1 4 


1 


-1 


41/120 


-1/20 


1/315 


-1/11880 


1/1330560 


6 


1 


-1 


1/4 


-1/42 


1/1008 


-1/55440 


1/8648640 



With Co = 0.973350, Ci = 0.005295, C2 = -0.006984, definitions ^ and ^ one can 
calculate partial structure functions again. The situation with the multipole L = 2 contribution 
is more complicated, this contribution appears to be rather significant especially for the large 
momentum transfer (the details can be found in ||48|]). The complete ^°^Pb partial structure 
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FIG. 12: Recoil energy dependence 
of the complete ■^''^Pb partial struc- 
ture functions from of 
Kosmas and Vergados [i^ . 
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functions Sfj^'^ from (|22j) of Kosmas and Vergados j48(| is given in Fig. as a function of the 
recoil energy. 

Concluding the discussion of lead structure functions we note, following |48(, that for a heavy 
nucleus and a high WIMP mass the momentum transfer dependence of the spin monopole 
(L = 0) matrix elements is quite large. It is, however, to a large extent neutralized by the 
spin quadrupole (L = 2). Thus the overall effect is not dramatic for WIMP masses less than 
100 GeV. For the spin-induced cross section of heavy odd-A nuclear targets, as in the case of 
^°''Pb, the reduction is less pronounced since the high multipoles tend to enhance the cross 
section as the momentum transfer increases (for LSP mass < 200 GeV) and partially cancel 
the momentum reduction 48^ . 



III. ON DATA ANALYSIS IN FINITE MOMENTUM TRANSFER FRAMEWORK 



To perform the data analysis in the finite momentum transfer approximation directly in 
terms of the effective spin nucleon couplings ao,i together with the scalar WIMP-proton cross 
section crgj(O) we propose to use formulas schematically given below. Nowadays it is also rather 
reasonable to assume crsi(O) ~ ^si(O)- '^^^ differential event rate dT)) can be presented as follows: 
dR{e, e) 



(25) 



T 



cpx Ma 
2m^ fil 



■{-)I{E^)d{E^-e)d{e-E^) 



(O'max) ^ 

[a'fUe^)]- 

4 / 

S'oo(g) + uj"^ 5*11 (g) + u Sm{q) 



2J+ 1 

°° (f 2) dv 



fiv)eiAfi\v'-2MAEn). 



JO (v) V 

Here (^max) = 4yU^ (i;^) where (v) and {v"^) are the mean and the mean squared velocity of 
WIMPs. The isovector-to-isoscalar nucleon couplings ratio is = ai/ao- Iii expressions (I25j) 
are introduced the detector threshold recoil energy e and the maximal available recoil energy 
e (e < -Er < ^)- In practice, for example with an ionization or scintillation signal, one has to 
take into account the quenching of the recoil energy, when visible recoil energy is smaller then 
the real recoil energy transmitted by the WIMP to the target nucleus. 

Formulas (j25|l allow experimental recoil spectra to be directly described in terms of only three 
independent parameters (cxgi, Og and u) for any fixed WIMP mass (and any neutralino 
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composition). Contrary to some other possibilities (see, for example, j43|), this procedure is 
direct and uses as much as possible the results of the most accurate nuclear spin structure 
calculations. 



IV. CONCLUSION 



There is continuous theoretical and experimental interest in existence of dark matter of the 
Universe. One of the best motivated non-baryonic dark matter candidates is the neutralino, 
the lightest supersymmetric particle. The motivation for supersymmetry arises naturally in 
modern theories of particle physics. In this work we discussed the spin-dependent interaction 
of neutralinos with odd-A nuclei. The nuclear structure plays an important role in determining 
the strength of the neutralino-nucleus cross section for this type of interaction and therefore 
defines the sensitivity of dark matter detection. In the limit of zero momentum transfer the 
relevant physical quantities are the proton and neutron spin averages (Sp(„)), which have to be 
evaluated within a proper nuclear model. These values determine the event rate expected in 
a direct dark matter search experiment due to spin-dependent neutralino-nucleus interaction. 
In our previous paper j?! the calculations of zero-momentum transfer spin-dependent matrix 
elements was reviewed and to our knowledge, a complete list of calculated spin matrix elements 
was presented for nuclei throughout the periodic table js^. The general feature is that spin 
matrix elements depend in general rather sensitively on the details of the nuclear structure. 
For a rather heavy WIMP (or light supersymmetric particle) and sufficiently heavy nuclei, the 
dependence of the nuclear matrix elements on the momentum transfer cannot be ignored. This 
affects the spin matrix elements. 

A comprehensive collection, to our knowledge, of the spin-dependent structure functions 
S^{q) calculated for finite momentum transfer (g > 0) within different nuclear models is pre- 
sented and discussed. These functions describe recoil energy dependence of the differential 
event rate due to spin-dependent neutralino-nucleon interaction, provided neutralino is a dark 
matter particle. Together with our previous paper "Nuclear spin structure in dark matter 
search: The zero momentum transfer limit" |3] this paper completes our review of the nuclear 
spin structure calculations involved in the problem of direct dark matter search. 

Now that spin structure functions are available for almost all experimentally interesting 
nuclei (and collected in this review), they could be coherently used by all experimental groups. 
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This will make easier and more reliable comparisons between results of different dark matter 
search experiments and put it on the equal footing I41)]. It will also allow one to reduce 
significantly the nuclear physics systematic uncertainties in the analysis of the data. 

This work was supported in part by the VEGA Grant Agency of the Slovac Republic under 
contract No. 1/0249/03 and by the Russian Foundation for Basic Research (grant 06-02- 
04003). V.A.B. thanks V.A. Kuzmin (JINR) for careful reading of the manuscript and his 
comments concerning the maximal momentum transfer. 



APPENDIX 

For completeness, in this appendix we collect some formulas which allow one to connect the 
WIMP-nucleon scattering with the WIMP-nuclear scattering. We directly follow Engel et al. 

Is^. The low-energy effective WIMP-nucleon Lagrangian is £efT = Xl^l^X ' <^m(^)' "where 
J^{x) Qc N'jfj^'-fr^N is the nucleon current. The one-nucleon matrix element of the current at 
finite q takes the approximate form 

{p,s\Mx)W, s') = U^ip, s) f ^^^^7.75 + ^¥^^.75) Ur^ip', s')e^'"'r (26) 

Here q^i = p^ — p'^-, Un{p, s) is the nucleon 4-component spinor and the energy transfer go was 
assumed to be very small. In the nonrelativistic limit the time component of current is 
proportional to v/c ~ 10"'^ and can be safely neglected. For the spatial component of the 
current one has the expression 

{p. s\M.)\p'. s') = (,|^i±^? - ilj^^^V'.-- (27) 

where \s) and \s') are two-component spinors. To obtain the cross section for scattering of 
the WIMP from nuclei one must evaluate the matrix element of the nucleon current between 
many-nucleon states. In the impulse approximation the cross section is 

(^;,g2) = I" I M = {s\a^\s') / d'x{J,M\J{x)\J,M') X e'^-^ (28) 



where |A1p is summed over s, s'M, M'. Here J is the angular momentum of the ground state 
and the nuclear current J^{x) is given by the sum over all nucleons with current matrix elements 
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from (j27|) . Expanding the current in vector spherical harmonics one obtains the form given in 
© and ©: 



i^,rj= or^7, ^81.(9)1/= E mwri'' mN)\' + \{N\\£l{q)\\N)r). 

\ ' ' L odd 



The transverse electric T'^^^iq) and longitudinal C^iq) multipole proiections of the axial vector 

nnnn 

current operator are given by (see [1|, 1221 |4Q|, |4]| for details): 



ri''{q) 



[-VLML,L+iiqri) + VL + lML,L-i{qri 



_1 \^(^ I 



VL + lMi,L+i(gr;) + \^ML,L-i{qn] 



where MLy^qri) = j l' {qri)[YL' {fi)ai]'" , is the pion mass and a^^^^^ is the isoscalar (isovector) 
effective spin-dependent WIMP-nucleon coupling. For example, the above-mentioned matrix 
elements one can calculate within the harmonic oscillator approach on the basis of results given 
in m. 
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